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ABSTRACT 

In the past, propulsion system simulations used In flight simulators have been extremely 
simple. This resulted In a loss of simulation realism since significant engine and aircraft 
Interactions were neglected and Important Internal engine parameters were not computed. More 
detailed propulsion svstem simulations are needed to permit evaluations of modern aircraft 
propulsion systems In a simulated flight environment. 

A real time digital simulation technique has been developed w^'ich provides the capabilities 
needed to evaluate propulsion system performance and aircraft system Interaction on manned 
flight simulators. A parameter correlation technique is used with real and pseudo dynamics 
In a stable Integration convergence loop. The technique has been applied to a multivariable 
propulsion system for use In a piloted NASA flight simulator program. Cycle time is 2.0 ms 
on a Unlvac 1110 computer and 5.7 n.s on the simulator computer, a Xerox Sigma 8. The model 
Is stable and accurate with time steps up to 50 ms. The program evaluated the simulation 
technique and the propulsion system digital control. The simulation technique and model used 
in that program are described and results from the simulation are presented. 

INTRODUCTION 

Cost Is a major factor In planning experimental ground and flight test programs. As a 
result, simulations, with their Inherent flexibility, are being used to a greater degree to 
design and analyze aircraft and propulsion systems controls before hardware is built and 
tested. 

Manned flight simulators have been used to evaluate aircraft stability and engine-out 
performance for various aircraft (1). The propulsion system models used in these simulation 
studies were extremi ly simple, however, providing only a thrust signal. This resulted In a 
loss of simulation realism to the extent that significant engine and aircraft interactions 
were not possible and important internal engine parameters were not available for analysis. 

lo overcome these deficiencies, reasonably detailed real time propulsion simulations are 
needed. Such simulations will provide the capability to evaluate propulsion systems and 
their Interaction with aircraft controls on manned flight simulators. The goal of the 
simulation development was to derive sucli a model and evaluate it in a piloted simulation 
program (2). This paper describes the development of a real time propulsion system 
simulation tliat provides this capability. The modeling technique was developed using the 
Under-the-Wlng (I'TW) version of the Quiet Clean Short-Haul Experimental Engine (QCREF.) and 
evaluated In the STOL aircraft de.scrlbed In (1). The simulation technique and propulsion 
system model are described and limited results from an evaluation of the simulation are 
presented and compared with experimental data. 

A non-linear propulsion system simulation such as (3) produces a model of high frequency 
fidelity running slower than real time In an analog format. Extension of this procedural 
development to hybrid models such as (4) produces a model with reduced high frequency 
fidelity but capable of Just real time. In a digital format as required for piloted 
simulators these models would require high sampling rates (small time steps) to maintain 
calculatlonal stability. Real time would be virtually impossible. The general approach 
taken In the real time digital simulator model development presented here was similar. For 
the level ol steady state and dynamic complexity required to meet this program objective, 
steady state accuracy does not have to be compromised over detailed models; but, high 
frequency content must be reduced significantly. The initial Improvement In execution time 


comes from efficient programming to attain minimum calculation time. More rapid Improvement 
occurs by maximizing stability to permit long time steps. 

R1:AL TIMK m.OTKD SIMUIJITION REQU1K1>1ENTS 

To satisfy the requirements of real time piloted simulation, Innovative matbema'ical modeling 
Is required. One must attain the desired level of fidelity yet have the computations 
accomplished In a limited amount of time. Also, since the propulsion system Is only a part 
of a larger simulation, only a fraction of the total computation time Is available for the 
pr Tulslon system calculations. Real time, then. In the context of overall simulation 
r. ,ulrements. Implies tliat the propulsion simulation must bo faster than real time to be 
• fectlve. 

Model Ing 

The basic technique used In fllgh'- simulation is to derive mathematical models that define 
aircraft system characteristics to the degree necessary to accomplish the objectives. The 
state of the aircraft system Is described by the equations of motion In terms of the aircraft 
accelerations, velc'citles and positions. Usually six degrees of motion are considered. The 
resulting differential equations are then Integrated with respect to time to obtain the 
aircraft states. Resulting kinematic Information Is then used as input to other simulator 
sub-systems such as visual motion, force-feel and instrumentation. The simulation equations 
are usually verv complex and are not amenable to analytical solution. Therefore the models 
must be programmed for numerical solution. Since the aircraft states can change rapidly with 
time, the computations are made often with small time Increments. A fast, accurate and 
stable Integration algorlthi Is essential. 

Real Time Monitoring 

The computer cannot meet real time requirements simply by repeating the calculation of the 
equations as rapidly as possible. This Is because the calculation time may be different for 
each time step. The calculations must be scheduled by a real time monitor to occur within a 
fixed time Interval that Is longer than the maximum time needed to solve the equations. Each 
Interval, oj frame time. Is measured by a real time clock In the computer. The computer Is 
interrupted at the end of each Interval and made to recycle through the equations. If the 
calculations cannot he co..ipleted In the .specified Interval a missed Interval will occur and 
real time operation may not be possible. At the end of each interval, time is Incremented by 
the frame time. This process, as shown at the top of figure 1, satisfies the real time 
calculat ion requirement . 

The real time requirement Introduces a conflict in objectives when selecting the Interval 
length. On the one hand the Interval must be minimized to retain higher frequency 
Information in the model. On the other hand, the model's detail and function are Improved by 
adding more equations which in turn incre.ises the calculation time. The conflict can be 
resolved by using multi-rate scheduling. 

Mult Irate Scheduling 

In multirate scheduling, equations are separated Into loops according to their frequency 
content as shown In figure I, Each loop Is run at a multiple of the basic frame time. In 
the case of three loops, the fast loop (high frequency) might be calculated four times as 
often as the slow loop (low frequency) and the Intermediate loop (middle frequency) twice as 
often as the slciw loop. The ratios can be vari.ihle. 

Wtille multirate scheduling permits higher frequency content In simulations, Innovative 
modeling techniques are still necessary to reduce computation times to a minimum and maintain 
simulation stability. The following paragraphs describe these modeling techniques and their 
application to a state-of-the-art turhofan propulsion system. 

rKOI’l'l.SU’.N SYSTf^l AITl ICATION 

The modeling technique was developed usli.g the I'nder-the-Wing (UTW) version of the Quiet 
(lean Short-Haul Experimental Engine (QCSEE) developed for NASA-Lewls by the General Electric 
I'ompanv (5) and shown In figure 2. 

The engine Is basically an EiOl core gas generator with a high bypass fan. It Includes a 
high M,ich Inlet, a variable pitch fan, a variable geometry fan duct exhaust nozzle and a 
digital electronic control system combined with a hydromechanical fuel control. The fan Is a 
li'w piessure ratio, low tip speed configuration with variable pitch blades and is driven by a 
low pressure turbine through reduction gears. 


The fan la capable of blade pitch changes from forward to reverse thrust. The fan pitch 
actuation and control are designed to move the blades from forward to reverse position In 
less than one second. 

The fan exhaust nozzle Is a hydraulically actuated variable area design. It is capable of 
area change from takeoff to cruise as well as opening to a flare position to form an inlet in 
the reverse thrust mode. 

The control system manipulates four variables to achieve rapid thrust response and noise 
suppression. Control of engine pressure ratio, fan speed and inlet mach number Is 
accomplished by manipulating fuel flow, fan blade pitch and exhaust nozzle area. Variable 
stator /anes are r’cheduled by core speed to attain optimum stall margin. 

The digital computer controls the output variables In response to consnands from the aircraft. 
It generates all control laws and logic and most of the limiting functions as well as power 
management, condition monitoring and failure indication and corrective action. The 
hydromechanical control provides an electrohydraullc servo fuel valve which Is used by the 
digital control for primary fuel control. It also provides backup fuel control through a 
core speed controller, acceleration and deceleration limits and primary control of the core 
compressor stators. 

modeling TECHNIQUE 

The analytical model Is derived from the real propulsion system. It represents 
mathematically the fundamental steady state and dynamic relations that exist between the 
engine components and controls. Engine dynamics are based on the dynamic form of the 
conservation equations and engine transient experience. Steady state performance is based on 
component representations derived from engine data. The control model is based on the 
control system specification. 

The form of the engine model and its Information flow are shorn schematically in figure 3. 
All major engine components are represented. The level of component detail is approximately 
the same as for state-of-the-art models. 

Since control evaluation is a prime consideration in propulsion system evaluations, a 
detailed control representation Is essential especially In the control law and switching 
logic areas. The control model is shown in figure 4. Accuracy of the control modal was 
a.ssured by deriving It directly from the control specification diagrams used In the real 
control. 

The dynamic engine simulation is derived bv . 1 ; plying ch*i basic conservation equations of 
continuity, energy and momentum to each component using the steady state compc'nent 
characteristics to define the boundary energy and mass across each boundary. Low frequency 
dynamics such as rotor speeds and component heat soaks are retained. High frequency dynamics 
such as volume dynamics are omitted. Algebraic loops that occur from omitted high frequency 
dynamics are converged using high gain integrators. The resulting dynamic effect Is similar 
to inductive lags. Component maps are generated by correlating input-output parameters to 
reduce complexity while retaining accuracy. These correlations are then curve-fit using 
segmented polynomial and geometric functions. The resulting analytical model is a twelfth 
order model which Includes four engine states, tour control sensor states and four control 
states. The resulting set of differential equations are solved as a general initial and 
boundary value problem using a two-step integration algorithm. The function generation and 
convergence techniques are innovative methods to accomplish rapid computation and stability. 

fTinction Generation 

Engine component performance maps are typically two or three variable functions. A fast, 
accurate function generation technique Is necessary to ensure real-time simulation. In this 
simulation technique, two approaches are used to achieve fast function generation. The first 
Is to simplify the basic functional relations through parameter correlation and the second is 
to approximate the functions as polynomial or geometric equations. 

This process will be described for the QCSEE variable-pitch fan map. The fan map, as 
determined from model data, Is a function of blade pitch angle, corrected speed and pressure 
ratio. This function was represented as shown in figure 5. The three-function map was 
approximated as a two variable basic fan map with multipliers on pressure ratio and corrected 
flow. The multipliers were functions of tlie blade pitch angle. The one and two variable 
functions were represented as multi-segment geometrical and polynomial functions. The 
functions were segmented to provide better accuracy over a wide range. The criteria used In 
developing each representation was to derive the relation with the least execution time 
consistent with a specified accuracy. The calculation time for the fan map function was 85 
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Ui conpartd to 7S0p ■ for a gancrAllzed function ganarator (6). Accuracy relallvc to the 
original fan map wae about one percent. Other multl-aegmant polynomial and analytical 
functions were calculated In about 15 to 30 pa. In thin way, all generalized function 
generation and subroutine calls were eliminated with substantial savlnga In calculation time. 

Convergence Technique 

Some form of convergence technique was required to solve the algebraic loop^ In the 
simulation. Pure Iteration was not used because of Its transient dependence. That Is, cycle 
or calculation time per time step varies with the rapidity of the transient. This is to be 
expected since non-steady mismatches occur during transients and more Iterations are required 
to converge the solution In order to avoid missed intervals. The frame time In piloted 
simulations which use pure Iteration must be based on the maximum calculation time even 
though the average calculation time may be considerably lower. Calculations! efficiency of 
Iteration Is therefore poor for real time simulation. 

In this simulation the iteration differences associated with algebraic loops were Integrated 
with high Integrator gains. A sketch of the process Is shown in figure 6. The procedure Is 
similar to using high gain integrators In analog computation to prevent algebraic loops. The 
previous "guessed" value, PCS, Is used to calculate the new required steady state value. The 
difference between these two values Is Integrated to generate an updated "guessed" value. 
The process shares the same Integration algorithm used to determine the model states. 

Referlng to figure 6, the Integration convergence method states that convergence takes place 
through a simple lag. In terms of a transfer function this Is:: 

PGS/P - l/(s/k + 1) 

where 1/K ” t. Tls the system time constant. The SMller the time constant (higher K) the 
faster the system response and the model approaches steady state convergence more rapidly. 
From a model convergence viewpoint then, high Integrator gains are deslreable consnensurate 
with stability requirements.. 

IWDEL EVALUATION 

The model wss exercised over a range of inputs to evaluate model stability and accuracy. A 
number of engine transients that might be expected to occur in a simulated flight e,raluation 
were run. 

Initially, the simulation was run to determine suitable Integrator gains and the effects of 
frame time on the results. A standard transient was used consisting of: (1) steady state 
operation at 62.5 percent, (2) a 62.5 to 100 percent power increase, (3) an in-to-reverse 
transient at 100 percent power and (4) an out-of-reverse transient at 100 percent power. 
Traces of net thrust and other variables were made to detect transient differences obtained 
with various combinations of gain and frame time. The traces were compared to those made at 
a gain of 100 and a frame time of 5 ms. This was done to insure that the data was compared 
to an accurate dynamic model. Selected results are shown in figure 7. 

Transient anomalies were expected with frame times larger than 50 ms since time constants In 
the engine and control models are higher than 50 ms. Time Constanta below 50 ms were 
purposely avoided In the model development. Deterioration beyond 50 ms was probably due to 
the single pass calculation of time constants In the control. Based on these observations 
the model was Judged to have a useful dynamic capability with frame times less than 50 ms. 

Every effort was made to eliminate unecessary calculations In the program. Divide 
computations were mlnlmlted and exponentiation eliminated. The model was programmed on a 
Unlvac 1100-40 computer using Fortran V. On that computer the model consumed 2.0 ms per time 
step. A somewhat simpler reverse model decreased that time to l.B ms. Cycle time was not 
transient sensitive. On the flight simulator computer, a Xerox Slgsa 8, cycle time was 5.7 
as or about three times slower. 

Transient Performance 

A number of transients were run which could be expected to be encountered In a mannev. ght 
simulation program to demonstrate the model's capability. A number of simulated engine and 
control component failures were also progransaed Into the model to analyze failure modes and 
effects In the flight program (2). A limited number of these transients are presented here 
along with corresponding experimental engine data acquired from recent testing at NASA-Lewis. 

Figure 8 shows a typical approach transient from approach power at 62.5 percent to go-around 
power at 100 percent. Shown are net thrust, fan speed and fan pitch angle. These are only a 
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few of the variables available from the model. Superlapoaed on the awdal traces are the 
experimental engine transient data. Discrepancies between these data are attributed mainly 
to the differences In the fan representation as derived from steady state cycle data and the 
real experimental fan performance as opposed to model inaccuracies due to siodel 
simplification. An updated fan map representation derived from the experimental data would 
eliminate most of the differences. 

The model accurately predicts the fan overapeed on acceleration. The difference In fen pitch 
angle and thrust are due to differences In Che fan pitch-fan flow relationship. The control 
Is manipulating fan pitch to maintain a scheduled fan speed. At the fan power required to 
maintain this speed, a higher fan pitch angle is required which In turn results In a lower 
fan flow rate and thrust. The effect Is really caused from lower actual fan performance than 
was predicted In the cycle deck model. 

Figure 9 showa the transient response of a simulated aborted takaoff-to-reverse aequence. 
Again, differences are attributed mainly to the difference In actual fan reverse performance 
and Che reverse map generated from the cycle deck. Of significance here are Che predicted 
thrust peaking on reverse Initiation and the predicted stall transition period. The model, 
however, does not accurately predict the actual fan speed droop during the stall transition 
to reverse. Again, this Is not due to Inaccuracies from model simplification, but due to the 
concept used In modeling stalled fan power consumption. 

CONCLUDING REMARKS 

The real time dlgitel propulsion system simulation developed under this program has generated 
a valuable simulation technology and flight simulator experience. It has yielded a feasible 
real time digital simulation approach. It has provided a fast, accurate and stable model for 
piloted flight simulation and application to other analytical controls problems. The 
technique.s developed so far have been effective in providing an adequate level of detail to 
evaluate propulsion systems In a simulaf'd flight environment. 
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Fqurt 1 - OCSEE engine modei schematic and information floN diagram. 

















Figurt 4. • CCSCE control modrl Khcnutic. 
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Figure 8. - Forward go-around transient. 
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Figure 9. - Aborted takeoff-to-reverse transient 


